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ABSTRACT 

We study contributions from inhomogeneous (patchy) reionization to arcminute scale 
(1000 < t < 10,000) cosmic microwave background (CMB) anisotropics. We show that inhomo- 
geneities in the ionization fraction, rather than in the mean density, dominate both the temperature 
and the polarization power spectra. Depending on the ionization history and the clustering bias 
of the ionizing sources, we find that rms temperature fluctuations range from 2 /zK to 8 /iK and 
the corresponding values for polarization are over two orders of magnitude smaller. Reionization 
can significantly bias cosmological parameter estimates and degrade gravitational lensing potential 
reconstruction from temperature maps but not from polarization maps. We demonstrate that a simple 
modeling of the reionization temperature power spectrum may be sufficient to remove the parameter 
bias. The high-^ temperature power spectrum will contain some limited information about the sources 
of reionization. 
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1. INTRODUCTION 

Free electrons in the reionized (or reionizing) intergalac- 
tic medium (IGM) with peculiar velocities lead to new 
anisotropies as the CMB photons scatter off of them. Here 
we revisit predictions for the temperature and polarization 
power spectra of these new anisotropies at small angular 
scales in light of recent indications from observations by 
the Wilkinson Microwave Anisotropy Probe (WMAP) for 
an extended period of reionization (Kogut et al. 2003). 

The argument for an extended period of reionization is 
as follows. The WMAP has detected the correlation be- 
tween temperature and polarization on large angular scales 
(Kogut et al. 2003) that has an amplitude proportional to 
the total optical depth of CMB photons to Thomson scat- 
tering, t (Sunyaev & Zeldovich 1980; Zaldarriaga 1997; 
Kaplinghat et al. 2002). Modeling reionization with a sin- 
gle sharp transition at z r i, a multi-parameter fit to the 
WMAP data gives z„ = 17 ± 5 (Spergel et al. 2003). 
On the other hand, the evolution of quasar spectra from 
z = 6.3 and z = 6.4 to z = 6 shows a rapid decrease 
in the amount of neutral Hydrogen, indicating the end of 
reionization (Fan et al. 2003). A simple interpretation to 
explain these two very different datasets is that reioniza- 
tion started early, z r i ~ 20, but did not conclude until 
much later (z ~ 6). 

As detected by WMAP, free electrons in the IGM scat- 
ter CMB photons and produce large angular scale CMB 
polarization. The fluctuations in the electron density field 
lead to small angular scale signatures both in the tem- 
perature and polarization. In the case of temperature, the 



new small scale anisotropy contribution is generally known 
as the kinetic Sunyaev-Zel'dovich (kSZ) effect (Sunyaev 
& Zeldovich 1980) Inhomogeneities in the electron den- 
sity may be due to perturbations in the baryon density or 
ionization fraction. If the ionized fraction is taken to be 
homogeneous and baryon density perturbations are calcu- 
lated to first order then the effect is called the Ostriker- 
Vishniac (OV) effect (Ostriker & Vishniac 1986; Vishniac 
1987). Non-linear effects from density perturbations make 
additional contributions beyond the OV effect on small 
scales (Hu 2000; Gnedin & Jaffe 2001; Ma & Fry 2002). 
When only modulations in the ionized fraction are con- 
sidered then we refer to the resulting temperature power 
spectrum as the patchy power spectrum. We focus on the 
patchiness because we find that, for the models we con- 
sider, variations in the ionized fraction are much more im- 
portant than variations in the baryon density. For a review 
of the imprint of reionization on the CMB see Haiman & 
Knox (1999). 

This extended period of reionization may be one in which 
the ionization fraction is highly inhomogeneous, i.e. "patchy" . 
As discussed by Aghanim et al. (1996), Gruzinov & Hu 
(1998) and Knox et al. (1998) (hereafter K98) such patchy 
reionization can lead to much larger anisotropies than ho- 
mogeneous reionization (see also Valageas et al. 2001, Agha- 
nim et al. 2002). Aghanim et al. (1996) and Gruzinov & 
Hu (1998) modeled the patches as spatially uncorrelated. 
K98 considered a model of correlated patches, and showed 
that these correlations were very important because they 
can greatly increase the signal in the £ — 1000 to 3000 
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range. In this paper we will take into account the corre- 
lations between these reionized patches by assuming they 
are associated with ionizing sources residing in dark mat- 
ter halos. The correlations will then be proportional to 
the linear theory dark matter correlations, with the pro- 
portionality constant given by an average of the halo bias 
weighted by halo mass and an ionizing photon production 
efficiency factor, since these determine the size of the indi- 
vidual ionized patch around each dark halo. In our models, 
the kSZ power spectrum typically increases by an order of 
magnitude when the inhomogeneities in the ionized frac- 
tion are taken into account. 

Patchy reionization also boosts the creation of polar- 
ization anisotropies on small angular scales. K98 pointed 
out that this polarization signal, since it is sourced by 
the temperature quadrupole as opposed to the velocities 
which have grown by gravitational instability, would be 
much smaller than the temperature signal. This was sub- 
sequently confirmed by explicit calculation (Weller 1999; 
Hu 2000). Hu (2000) pointed out that the B-mode power 
spectrum would be equal to the E-mode power spectrum 
(for a flat-sky version, see Baumann et al. 2002). We also 
calculate the polarization power spectrum from reioniza- 
tion and find that it is indeed too small to significantly 
affect cosmological parameter determination or lensing po- 
tential reconstruction. 

Our paper is organized as follows. In § 2 we describe 
the kinetic SZ effect and how to calculate the temperature 
and polarization angular power spectra in the presence of 
ionization fraction and baryon density fluctuations. In § 3 
we briefly describe our models for reionization (described 
with more detail in Haiman & Holder 2003) , including our 
calculation of the effective bias of the patches. In § 4 we 
discuss the resulting power spectra and how they can be 
used as probes of reionization. We also derive in § 4.1 
the expected errors from a ground-based experiment ca- 
pable of detecting this signal. In § 5 we calculate the 
resulting bias in the cosmological parameter determina- 
tion from Planck data that can be caused by patchy reion- 
ization, and show how simple phenomcnological modeling 
can eliminate these biases. Finally, in § 6, we discuss the 
impact on lensing potential reconstruction. We conclude 
with a summary in § 7. Throughout the paper, unless 
stated otherwise, we will assume the cosmological param- 
eters, il m = 0.29, n A = 0.71, n b = 0.047, h = 0.72, 
t = 0.17 and n = 1 with normalization erg = 0.9 in agree- 
ment with a recent analysis of CMB data (Spergel et al. 
2003). 

2. KINETIC SZ EFFECT 

Compton-scattering via moving electrons leads to a Dop- 
pler effect and a corresponding temperature perturbation 
along the line-of-sight: 



AT f 
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where n is the unit vector along the direction of obser- 
vation, r\ is the conformal time (with value rjo today), 
x = n(?7o — fj) and v is the electron velocity. The visibility 
function, g(x,r]), corresponds to the probability that the 
photon scattered at time r\ and traveled freely since then: 

a T n e (x, 77) e~ T( - r,) , 



where r is the optical depth to electron scattering between 
here and distance r] — r\ and ot is the Thomson cross 
section. The free electron density is n e = X e (x, ry)n p (x, rf) 
and both the ionization fraction, X e and proton density, 
n p , are allowed to have a spatial dependence. 

Under a smooth distribution of free electrons with n e = 
n e , the Doppler contribution to temperature anisotropies 
is significantly reduced as photons scatter against the crests 
and troughs of the velocity perturbations (Kaiser 1984). 
This cancellation is avoided when there are significant fluc- 
tuations in the free electron density field, which can orig- 
inate from perturbations in the ionized fraction (6x e ) or 
the baryon density (5b). The total contribution to the 
temperature anisotropy, up to second order, can then be 
expressed as: 

AT f 

— {n)=a T n p {r l0 ) / dr] [gT 2 : (r7)e~ T ' X e (rjj\ nq, (3) 
with 

q = (1 + s Xe + s b ) v (4) 

and cr T n p (?7o) = 2.03 x 10~ 5 tt b h 2 Mpc" 1 . 

2.1. Power spectrum 

We next calculate the kSZ power spectrum on small 
angular scales, where the signal is not overwhelmed by 
the primary anisotropies. Then, for £ ^> 1, we can use 
the Limber approximation (Limber 1954) and write the 
spherical Fourier coefficient C\ in terms of the 3D power 
spectrum of q(x, 77) (Kaiser 1992; Dodelson & Jubas 1995; 
Hu & White 1996; Jaffe & Kamionkowski 1998): 

C e = cr^npivo) 2 / ^-a~ 4 (?7)e~ 



(5) 

where 

(OL(k)-ql(k')) = 2(2n) 3 5(k-k')P q Jk) (6) 
and q^ (k) is the component of the 3D Fourier transform 
of q(x) that is perpendicular to k (qj_(k) - k = 0). Thc_im- 
plicit assumption here is that the quantity [a~ 2 (r])e~ T X e (r])\ 
is slowly- varying across a wavelength of the perturbation 
(A?7 <~ i]o/£ at most). This is a good approximation for 
the angular scales of interest (£ <; 1000). 

The key quantity to compute is therefore P q± for all 
possible contributions. To first order q = v and v(k) oc k 
so that, as already stated, there will be no contribution 
to the CMB power spectrum from equation (5). At larger 
angular scales, on the other hand, the anisotropy from the 
component of q(k) parallel to k becomes significant and 
is taken into account in standard Boltzmann calculations. 
Considering the density-modulated and ionized-fraction- 
modulated anisotropies, the components of q(k) perpen- 
dicular to k are 

[&(|k-k'|) + fa,(|k-k'|)], 

where fi' = k • k' and we have assumed that the velocity 
field is a potential flow, v(k) = — iV(k)k, as expected 
in linear perturbation theory. The corresponding power 
spectrum is then 

1 f d 3 k' 
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Note that, due to geometric reasons, there is no connected 
fourth moment contribution to P q± (k) (Cooray 2001; Ma 
& Fry 2002). 

2.2. Density-modulated anisotropies 

If reionization was homogeneous throughout the Uni- 
verse, temperature anisotropies can only arise from pertur- 
bations in the baryon density (e.g. a = b = 5b in equation 
8). When 5b is in the linear regime, this is the Ostriker- 
Vishniac effect (Ostriker & Vishniac 1986; Vishniac 1987), 
which has already been studied in some detail (Dodclson 
& Jubas 1995; Hu et al. 1994; Persi et al. 1995; Hu & 
White 1996; Jaffe & Kamionkowski 1998; Hu 2000). Here 
we will just try to give an order of magnitude estimate for 
the effect. For wavelengths that are much smaller than the 
coherence scale of the velocity fields, equation (8) simpli- 
fies considerably. First we note that, using the continuity 
equation: 

P vv (k',ri) cx P 5c sAk')/k' 2 Ps c v{k',v) « P Sa s c (*')/*', 
where we have assumed that the baryonic gas traces the 
cold dark matter. For large k', typically larger than k m ~ 
0.1 ft/Mpc (Hu 2000) we then have P vv ~ and P ScV ~ 
(density fluctuations uncorrelated with the bulk velocity 
field). Therefore, for large k the integral in equation (8) 
can be evaluated under the approximation |k — k'| w k so 
that, 

1 



PqAk'V) ~ 7;Ps c 8 c (k,r))v rms (r)), 



with 
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The factor of 1/3 arises from angular averaging, since 
only the line of sight component of the velocity field con- 
tributes to the signal (the two transverse parts drop out). 
For these small wavelengths, typically Ps c s c (k) cx fc -3 so 
that, introducing the above expression in equation (5), 
gives d cx £~ 3 (for £ ^> i] k m ). Moreover, assuming our 
LCDM fiducial model and instantaneous reionization of 
the Universe, we obtain, as an order of magnitude esti- 
mate, e 2 C t /(2TT) - 3 x 10- 13 (5000/£). 

In the low redshift universe most of the baryonic mat- 
ter resides in collapsed, high density objects. We there- 
fore need to modify the above expressions to take into ac- 
count the non-linearities of the density perturbations (Hu 
2000; Gnedin & Jaffe 2001; Ma & Fry 2002). These non- 
linearities only affect the density field below the coherence 
scale of the bulk velocity. Therefore, equation (9) will be 
a good approximation for scales that become non-linear. 
One then only needs to obtain the corresponding non- 
linear matter power spectrum using an appropriate pre- 
scription such as the halo model (Ma & Fry 2000). Figure 
1 compares the OV power spectrum with the non-linear 
correction calculated using the halo model (thin lines). 
We see that for sufficiently large £ there is an increase in 
power due to a similar increase in the matter power spec- 
trum. We point out that this non-linear OV effect includes 
contributions to the kSZ signal from low redshift galaxy 
clusters. As will be seen below, the amplitude of the kSZ 
signal from the high-redshift ionized patches overwhelms 




Fig. 1. — The several contributions to the kSZ temperature power 
spectrum assuming the LCDM model described in the introduc- 
tion. The thin solid line is the O-V power spectrum from homo- 
geneous reionization and dashed line is the non-linear correction. 
The thick solid line is the typical contribution from inhomogeneous 
reionization and uses the model shown in Figure 3 and 4 (solid 
line). Dashed line shows the approximate model from eq. (17). 
The patchy power spectrum basically follows the shape of the mat- 
ter power spectrum. Also shown is the primary lensed CMB power 
spectrum (dot-dashed). 



this low-redshift contribution in most reionization mod- 
els. The results in Ma & Fry (2002) agree reasonably well 
with hydrodynamical simulations (Springel et al. 2001; da 
Silva et al. 2001; White et al. 2002) for the angular scales 
of interest. 

2.3. Patchy reionization 

Exactly as with perturbations in the matter density, 
inhomogeneities in the ionized fraction create tempera- 
ture anisotropies through modulation of the Doppler effect 
(K98, Gruzinov & Hu 1998). Considering the case where 
a = b = 5x c in equation (8), we now need to obtain an ex- 
pression for Pg x s x (k)- We assume the ionizing radiation 
is coming from stars formed from gas clouds that cooled 
in dark matter halos. The dark halo correlation function 
can be expressed in terms of the linear theory dark matter 
one (Mo & White 1996; Jing 1998): 

Z halo (r;M,M',r))=b(M,r])b(M',r))(6 aSc (r,r)), (U) 
where £,haio(r; M, M') is the correlation function between 
halos of mass M and M' and 6() is the clustering bias pa- 
rameter. To obtain the clustering properties of the ionized 
fraction we need to take into account the dependence of the 
size and spatial distribution of the ionizing patches on the 
masses and ionization efficiencies of dark halos. Therefore, 
one should be able to write the correlation as: 

where b e f / is the mean bias weighted by the different halo 
properties. Since these ionizing patches are rare objects 
we expect them to be more clustered than the underly- 
ing dark matter (e.g. 6 e // > !)■ We will compute this 
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bias in more detail in Section 3. Finally, one has to take 
into account that the perturbations will be smoothed out 
on scales smaller than the size of the patch. Choosing a 
Gaussian filter, the power spectrum for the ionized fraction 
will then be 

Ps x . Sx .(k, V ) = b 2 eff (r])P 5c5c (k,v)e- k2R \ (13) 
where R is the mean radius of the ionized patches (HII 
regions). 

We crudely model the time-dependence of R as 



R- 



1/3 



R. D 



(14) 



1-X e 

where R p is the (comoving) size of the fundamental patch 
which we take to be R p <~ 100 Kpc. The cutoff scale R 
is equal to R p when reionization starts and increases with 
time as HII regions overlap and form larger HII regions. 
Equation 14 follows from assuming a Poisson distribution 
of patches and that their volumes add when they overlap. 
But the essential consequence of equation (14) is that on 
the scales of interest (£ < 10, 000 or R > IMpc) the tran- 
sition is very sudden; i.e., R increases from 1 Mpc to oo 
very rapidly. We expect this behavior to be generic and 
thus our results should be insensitive to the specifics of the 
overlap modeling. 

To calculate the patchy power spectrum we neglect the 
cross-correlation between velocity and 6x e since it is small, 
though we discuss this contribution later. Then, we simply 
plug expression (13) into equation (8). The result for our 
fiducial reionization model (discussed in the next section) 
is shown in Figure 1. To have a better understanding of 
the expected shape and amplitude of this power spectrum 
we will now derive an approximate expression. Using an 
expression similar to equation (9) we find the patchy power 
spectrum to be 
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where r\i signals the start of reionization (when X e starts 
increasing from its value at recombination) and rj\ corre- 
sponds to the time when X e = 1. We can write the time 
dependence of the dark matter perturbation as S c (k, -q) = 

(<5o ) ^ c (^' %)> wncre G is the growth factor (Carroll ct al. 
1992) and Go is the value today. Moreover, assuming 
Xi — x\ <C xi, we can pull the I dependence out of the 
time integration, so that: 



where G 



-Pn 

f 



drj a 



Go Go 
X 2 e (rj)bt ff (ri)(GG) 



(16) 



dG 
drj 



The quantities, x, R and r are assumed 

to be evaluated at the average time, (771 + rji)/2. At the 
rcdshifts of interest for the above time integration, the 
cosmological constant is negligible, so that G ~ a and H <~ 
Hoy^flm/a 3 . The power spectrum can then be expressed 
as 



6.18 x 10 

dz 



2l eT 2 (l/x)e 

-2/ x ,2 



-.Xi(z)K /f (z), 



(17) 



where we are using our fiducial cosmological model, so 
that v rms (r) )/c = 1.8 x 10~ 3 , x ~ 8.0 x 10 3 Mpc/h and 
R ~ 0.20 Mpc/h. The matter power spectrum is nor- 
malized to give as = 0.9 and the transfer function, T(k) 
(Bardeen et al. 1986; Eisenstein & Hu 1998; Seljak & Zal- 
darriaga 1996), is normalized to 1 on large scales. All 
the dependence on the reionization model is in the time 
integration, through X e and 6 e //. Therefore we expect 
the power spectrum from different reionization histories 
to have approximately the same shape but different am- 
plitudes. Figure 1 shows a possible power spectrum from 
the above patchy reionization model. Also shown is the ap- 
proximate model from equation (17). This signal typically 
overwhelms the OV effect, even when including the non- 
linear correction, which can be traced back to the large 
effective bias, 6 e // <; 5, in equation (16). 

We have ignored the term, (wi<5x e2 )( w 2'5x cl ) as sub- 
dominant to (5x e i5x e 2)( v ^ v ^}- The former term may be- 
come important for length scales above the velocity coher- 
ence length, which projects to I ~ 500 from high rcdshifts. 
Due to the tendency of overdense regions to fall toward 
each other, this term is negative and the net effect will be 
a suppression of power. This power suppression will be 
small in the scales of interest. We have also neglected the 
cross correlation (6 c 6x e ) because it will be suppressed by 
approximately one factor of the bias from (Sx c SxJ- 

2.4. Polarization Signatures 

Thomson scattering of radiation with a quadrupole ani- 
sotropy generates linear polarization in the CMB. Again, 
the dominant source at small angles is the modulation of 
the quadrupole through perturbations in the matter den- 
sity or ionized fraction. The power spectra of the E and 
B parity states are then (Hu 2000; Baumann et al. 2002) 



r>w _ riot, _ 



3c4 n p (r) ) 2 f drj 



100 J :>~ 

l2 



■a-\rj)e-*Xi{rj) x 



{P6 c 6 c +P6 Xe 5 Xe ) Q 2 rms (v), (18) 
where Q r ms is the quadrupole anisotropy. The power spec- 
trum from the primordial quadrupole source will typically 
dominate over the kinematic and intrinsic contributions 
(Hu 2000). Moreover, this quadrupole is a slowly varying 
function of time and can be taken out of the integral above. 
Figure 2 shows the patchy and density contributions to po- 
larization (using Qrms ~ 25/iK), compared to the B mode 
polarization from tensor and lensing effects. Again the 
patchy signal dominates over the density modulated po- 
larization. However, it does not seem to be an important 
contaminant on the relatively large scales (I <> 200) where 
tensor B modes are significant, and it is well below the 
scalar E mode on all scales. 

The patchy polarization power spectrum is <~ 10 5 times 
smaller than the patchy temperature power spectrum be- 
cause Q rms ~ 10 _2 w rms and due to the geometric factor of 
100 in equation 18 in place of the factor of 3 in equation 15. 

3. REIONIZATION MODELS 

We now_ need to obtain the evolution of the electron 
fraction, X e and effective bias, b e ff, in physically moti- 
vated models of the reionization history. For that purpose, 
we adapt semi-analytical models from Haiman & Holder 
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Fig. 2. — B mode polarization for our fiducial model, with contri- 
butions from lensing (thin solid line) and tensor modes (thin-dashed) 
with a tensor to scalar ratio of 10 - . Also shown is the contribution 
from density (thick dashed) and ionization (thick solid) modulated 
scattering. The density modulated contribution assumes gas traces 
dark matter and uses the halo model for non-linear corrections. Also 
shown for comparison is the scalar E mode contribution (dot-dashed 
line). 



ements from the first generation. This population contin- 
ues ionizing hydrogen, but is also self-limiting: gas infall 
to these relative shallow potential wells is prohibited inside 
regions that had already been ionized and photo-heated 
to 10 4 K (Thoul & Weinberg 1996). As a result, hydrogen 
reionization starts slowing down around z <~ 10. However, 
at this stage, still larger halos with virial temperatures of 
T^2x 10 5 K start forming. These relatively massive halos 
are impervious to photoionization feedback, and complete 
the reionization of hydrogen at around z w 7. 

In order to predict the kSZ signal, we need to compute 
in our models the mean bias of the ionized regions as a 
function of redshift. At any given redshift, there is an 
ensemble of HII regions, the sum of which make up the 
total ionized volume fraction. Each HII region has at least 
one ionizing source, each of which has a fixed halo mass 
and formation redshift. 

We assume that the volume of each ionized region scales 
linearly with halo mass within each of three distinct mass 
bins; we only allow a dependence of ionizing photon pro- 
duction efficiencies in discrete jumps between these bins. 
This assumption could be relaxed to include more generic 
dependencies on halo mass. This, in general, would change 
the effective bias we obtain, but would involve introduc- 
ing more free parameters; we avoid this complication in 
the present paper. Under these assumptions, the mean 
bias of ionized regions can be computed as follows. We 
first note that at a given redshift, the total volume-filling 
fraction Frit is given by 



(2003), modified to include the spatial biasing of the ioniz- 
ing sources and hence of the ionized regions. These mod- 
els fit all the relevant observations (including the electron 
scattering optical depth r = 0.17 measured by WMAP 
and inferences from redshift z ~ 6 quasar spectra). For 
details, the reader is referred to that paper; here we pro- 
vide only a brief summary, and a description of how the 
bias is calculated. 

We follow the volume filling fractions .Fhii of HII regions 
(which, by definition, equals the mean electron fraction 
X e ) assuming that discrete ionized Stromgren spheres are 
being driven into the IGM by ionizing sources located in 
dark matter halos. In this picture, each fluid element is 
engulfed by an ionization front at a different time. Reion- 
ization has a complex history that reflects contributions 
from three distinct types of ionizing sources, and two dif- 
ferent feedback effects (all of which have physical motiva- 
tions as described in detail in Haiman & Holder 2003). In 
short, ionizing sources (assumed to be massive metal-free 
stars) first appear inside gas that cools via H 2 lines, and 
collects in the earliest non-linear halos with virial temper- 
ature of 100 K <, T & 10 4 K (Haiman et al. 1996). These 
sources ionize ~ 50% of the volume in hydrogen. How- 
ever, at this stage (redshift z ~ 17) the entire popula- 
tion of these sources effectively shuts off due to global H 2 - 
photodissociation by the cosmic soft UV background they 
had built up (Haiman ct al. 1997a, b, 2000). Soon more 
massive halos, with virial temperatures of 10 4 K & T ^ 2x 
10 5 K form, which do not rely on H 2 to cool their gas (they 
cool via neutral H excitations), and new ionizing sources 
turn on in these halos. These are assumed to be "normal" 
stars, since the gas had already been enriched by heavy el- 



-Fhii (z) = Pb(z) 

rf-Fcoll.Ia 



dz i eib ; (z ) + [I 



dz 



dz 



, i \ . dF co \\ u , 
{z)+eil ^z- {z) 



Fmi{z')] x 
Vun(z',z), (19) 



where /?b = ^bPcrit is the average baryonic density and 
p cr it is the critical density of the universe. In the sec- 
ond term on the right hand side of equation (19), we have 
explicitly included a factor (1 — Fun), which takes into ac- 
count the photoionization feedback. Here i^coiiji, ^coii.ia 
and Fcoii,ib, are the fractions of baryons collapsed into ha- 
los of increasing mass ranges. For example (using p to t as 
the present-day total mass density), 

i /-M max . 

Foiui = — / dM^-{z)M (20) 

where M m ; n and M max correspond to halo virial tempera- 
tures of XUr = 100K and T v ; r = 10 4 K at redshift z, respec- 
tively. The en, ei a and eib are the efficiencies in these halos 
of producing and leaking ionizing photons into the IGM. 
The efficiencies are defined as the product e* = N^f*/^, 
where /* = M Jf /{Q,\ ) M\ lg \ /Q, m ) is the fraction of baryons 
in the halo that turns into stars; A 7 is the mean number 
of ionizing photons produced by an atom cycled through 
stars, averaged over the initial mass function (IMF) of the 
stars; and / esc is the fraction of these ionizing photons that 
escapes into the IGM. Finally, Vnu(z n, z) is the volume 
per unit mass and unit efficiency ionized by redshift z by 
a single source that turned on the earlier redshift z OD > z. 
The evolution of Fun for a few models described in Haiman 
& Holder (2003) is shown in Figure 3. They all produce 
a total electron scattering optical depth of r = 0.17 by 
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Fig. 3. — X e for three different fiducial models with r = 0.17. 
Solid line is for the fiducial model used throughout the paper. 
Dashed line assumes that H2 feedback effectively shuts off star for- 
mation in all Type II halos. Dot-dashed line assumes that a transi- 
tion from metal free to normal star formation takes place abruptly 
at a certain redshift. 



construction, in agreement with the WMAP result. For 
the one we take as our fiducial model in this paper (solid 
line), the efficiencies are en = 200 and ei a = eib = 80. 

To calculate the effective clustering bias for any field, 
the first step is to write the field as an integral over halo 
mass and formation rcdshifts, which we have effectively 
done with equation (19) and (20). We then use this sum 
to define 6 e // as the appropriate weighted average of the 
bias factor for a halo of mass M and formation redshift z, 
b(z,M). The effective bias for the ionization fraction for 
our reionization models can thus be efficiently written as: 

b eff (z) = F& u /F ml ( 21 ) 

where 

F&n (z)^p h (z) J Z dz'^e lh b 1 ^)^^(z') + 

/ /\ ^-^coll.Ia / /\ , i / l\ ^-^coll.II / i\ 

ei a 0la(2 ) ; — ■ — \Z ) + £11011(2: ) ; (Z ) X 

az az 

[l-F mi (z')]^V mi (z',z) (22) 

and 

b n {z) = [Fcoii.iiPtot] / dM—(z)b(z,M)M, 

with analogous expressions defining b\ a {z) and bih(z). 

In this last equation, b(z, M) is the usual halo bias (Mo 
& White 1996). The averaging over formation redshift 
only introduces a small correction to the bias calculation, 
since at a given redshift most of the contribution comes 
from young ionizing sources. 

Figure 4 shows the effective bias for the same three mod- 
els as in Figure 3. The increase of the bias as z decreases 



Fig. 4. — Effective bias (equation 21) for the same fiducial models 
as in Figure 3. 

toward ~ 7 is due to the formation of rare, massive halos, 
which play an important role in the last stages of reion- 
ization. 

Note that the bias factor as calculated in Eq. (21), Eq. (22) 
and (23) agrees with the halo model description of source 
bias (Cooray & Sheth 2002) when one makes the connec- 
tion between the presence of ionized electrons and the halo 
occupation distribution. Here, the bias factors are calcu- 
lated with the assumption that ionized electrons occupy 
halos with mass in the range M m - m and M max , such that 
the mean halo occupation number is unity in this mass 
range and zero otherwise. 

We now remark on the robustness of the model of patchy 
reionization we have presented. This simple picture has 
the HII regions localized around the halos hosting the ion- 
izing sources. Further, our calculation of V(z, z 1 ) assumes 
the recombination rate is given by the square of the back- 
ground density times some constant clumping factor. In 
reality, there must be complicated shapes for the volume 
V(z,z'), which depend on the three-dimensional arrange- 
ment of the density field, and also on the locations and 
strengths of the ionizing sources. The degree of localiza- 
tion of the HII regions around the sources of ionizing pho- 
tons will also depend on how hard the radiation is, with 
harder spectra reionizing more uniformly (Oh et al. 2001; 
Venkatesan et al. 2001). And, if reionization is from de- 
caying dark matter then X e would be yet more uniform 
(Hansen & Haiman 2003). 

Another simple, but complementary, picture (Miralda- 
Escude et al. 2000) (hereafter MHR) is that the ionizing 
radiation from a source leaves the surrounding dense re- 
gions neutral and the i-fronts propagate along directions of 
the the lowest column density into the voids. In this case 
it actually does not matter that the sources are in high-cr 
peaks. MHR postulate that all gas below some threshold 
density, p tl is ionized, and all gas above this pt remains 
neutral. The patchy power spectrum in this case could 
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Fig. 5. — Patchy power spectra for the reionization models in 
Figure 3 (same line styles), together with other astrophysical con- 
tributions and expected measurement errors (see text). The solid 
(dashed) straight line is the the point source contribution at 217 
GHz before (after) multi-frequency cleaning. The primary unlcnsed 
(dashed) and lensed (solid) CMB power spectra are also shown as is 
the thermal SZ power spectrum from White et al. (2002) (dotted) 
with its expected amplitude at lower frequencies. The thin line close 
to the solid one shows the patchy power spectrum for a model with 
T = 0.11 but large bias. 

be computed using the two-point function of regions with 
P < Pt- 

The MHR picture is the extreme opposite of ours (see 
also Peebles & Juszkiewicz 1998). Which picture is a bet- 
ter description of reality depends entirely on how rapidly 
percolation occurs. We assume it only happens as X e gets 
very close to unity. For our calculations to be correct at 
I < 10, 000 requires that prior to X e ~ 1, the ionizing radi- 
ation in the HII regions mostly comes from sources closer 
than 3 Mpc. 

4. TEMPERATURE POWER SPECTRA 

In Figure 5, we plot the patchy power spectrum for the 
three reionization models (including our fiducial model), 
all with the same cosmological parameters and r = 0.17 
but differing X e (z) and b e g(z). As expected, the shapes 
are very similar. Allowing r to vary also does not sig- 
nificantly affect the shape, as demonstrated by the fourth 
model plotted with r = 0.11 but higher bias that is close to 
a r = 0.17 model. The higher bias is achieved by moving 
the ionizing sources to more massive halos (Type la and 
lb). The optical depth is prevented from dropping even 
further by increasing the efficiencies in Type la halos. 

The degeneracy between t and bias means that r can- 
not be inferred from the kSZ power spectrum alone, as 
suggested by Zhang et al. (2003). However, if r becomes 
tightly constrained by polarization measurements at low £, 
which is possible with future missions (Holder et al. 2003), 
the mean effective bias for the sources of reionization can 
be determined. Having this information will constrain the 
typical halo masses of the sources that dominate reioniza- 



FlG. 6. — Patchy power spectrum for variations in r about two of 
the models shown in Figure 5 (thick solid and dashed lines). Solid 
lines correspond to the following values of r (from bottom to top): 
0.05, 0.11, 0.17, 0.24, 0.31. Dashed lines correspond to r =0.11, 
0.17, 0.22, 0.28. Also shown are the expected measurement errors 
and primary lensed CMB power spectrum. 

tion. Whether the sources are located in the 2-sigma or 
in the 3-sigma peaks is important for the physics of reion- 
ization (see Rees 1999). This will also constrain the mean 
efficiency e. 

To demonstrate the range of possible power spectrum 
amplitudes we calculate the patchy power spectra for a 
variety of reionization models with values of r spanning 
the 2cr range allowed by WMAP. In Figure 6 we plot the 
patchy power spectrum with variations in r about two of 
the models shown in Figure 5 (solid and dashed lines). 
Variations in r for the solid (dashed) line are done by es- 
sentially changing the efficiencies in the Type II (la) halos. 
Again the shape is very similar (at least up to I ~ 10 4 ) 
with the amplitudes ranging from 1 x 10~ 12 to 8 x 10" 12 
at £ ~ 2000. 

At I > 10 4 the power spectrum is sensitive to the struc- 
ture of the reionized patches on comoving length scales 
smaller than 1 Mpc. At these smaller scales the details of 
the patch geometries and sizes become important. With 
sufficiently reliable modeling of the power spectrum on 
these scales (which we have not attempted) observations 
at I > 10 4 could be used to constrain these patch prop- 
erties, although contamination by point sources will make 
such observations difficult. 

4.1. Measuring the kSZ Power Spectrum 

Ground-based observations, with 1' to 2' angular resolu- 
tion, may be able to measure the kSZ power spectrum. In 
order to do so, several other astrophysical signals must be 
cleaned from the maps: namely, thermal SZ and emission 
from dusty galaxies. 

We take as an illustrative example a multi-frequency ex- 
periment that includes a channel at 217 GHz. We assume 
that the thermal SZ in the 217 GHz channel is insignif- 
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icant, or at least that the residuals after multi- frequency 
cleaning are insignificant. We further assume the 217 GHz 
channel has an angular resolution of 0& = 0'.9 and an rms 
error of At = 12/uK on every beam-size pixel in a fraction, 
/sky = 0.01, of the sky, or 400 sq. degrees. 

For our point source model, we adopt model E of Guider- 
doni et al. (1998). For this model, with 0'.9 resolution, 10 
mJy point sources stand out as 5cr fluctuations where a 2 
is the variance due to sources fainter than 10 mJy. Re- 
moving these bright sources with a 5c cut leaves a power 
spectrum due to the faint point sources of l 2 Cf s /(2n) = 
2.5 x 10~ 12 (//2000) 2 and is plotted as the solid straight line 
in Figure 5. Exploiting the frequency dependence of these 
remaining sources (with respect to CMB temperature fluc- 
tuations) might reduce the power in this contamination by 
a factor of 10 (dashed blue line). 

The error on the total fluctuation power is 



AC, tot = 



-(Cf* +N t ) where 



^ (2Z + l)/ sk / 

Cf * - Cf + Cf s + Cf sz and 
AT; = (A T 6 ) 2 exp(Z 2 2 /(81n2)). 



(24) 



Here Cf corresponds to the primary anisotropy contribu- 
tion and C ; kSZ is the power spectrum from the kSZ effect. 
A multi-parameter fit to C kSZ , the amplitude of a shot- 
noise spectrum, and the cosmological parameters govern- 
ing Cf will result in 



AC, 



■kSZ 



AC, 



(25) 



since Cf and Cf s will be very well determined. We plot 
the AC, tot errors in Figure 5 about our fiducial kSZ model. 

5. IMPACT ON PARAMETER DETERMINATION 

The large anisotropics at t <; 1000 predicted from a 
patchy reionization phase may significantly affect the high 
precision parameter estimation predicted for future CMB 
experiments. Let us suppose we obtain a set of best fit 
parameters a, by minimizing the likelihood, £° , without 
considering the patchy reionization signal (e. g. 91 g a C ° = 
0). The above set should be shifted by an amount Scii so 
that they minimize the true likelihood, C n , which includes 
the patchy contribution: 



6a, = J2 ^ij 



dan 



where Tn* is the fisher matrix, 



d 2 ln£ n 

dajdai 



(26) 



(27) 



The bias in the estimation of each parameter will then be, 
as in K98, 



8a; 



E iCf^dCf^ 



(28) 



it 



where Cf can represent the power spectrum from inhomo- 
geneous reionization and at is the error expected on the 
Ce's, e.g. the sum of the errors due to sample variance 
and noise, Ng (Knox 1995; Bond et al. 1997; Zaldarriaga 

1997), ^2m(C[ + Cf + N e ). We estimated the bias for 
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Fig. 7. — Absolute ratio of bias to la for Planck as a function of 
(max- Shown are the parameters that will have a significant bias: 
n s (continuous line), m v (dotted), Yu B (dashed), n' a (dot-dashed) 
and w x (long-dashed). 



the same set of cosmological parameters as in Kaplinghat 
et al. (2003) taking into account the weak lensing contri- 
bution on small scales. Using our fiducial model for Cf, 
Tabic 1 shows the expected bias for Planck (Tauber 2001). 

In Figure 7 we plot the ratio of this systematic offset to 
the statistical uncertainty as a function of the maximum 
multipolc considered. The situation is even worse for a 
sample variance dominated experiment. As we can see in 
Tabic 2, most of the biases are already important at I = 
2000, thus becoming a problem for parameter estimation. 

The solution is to obtain a model for the patchy power 
spectrum in terms of a few parameters and take them 
into account when doing a multi-dimensional fit to the 
data. Looking at equation (16) we see that although the 
amplitude depends on the reionization model, the shape 
does not. Let C\ be a typical Cg from patchy reioniza- 
tion, then any other reionization model can be written as 
C{ w AC}. There will still be a bias due to the slight 
differences in shape between the reionization models (e.g. 

Cf — Cl — AC}). Table 3 shows the ratio, bias/a, for a 
cosmic-variance dominated experiment when considering 
the amplitude as an extra parameter. We see that the ex- 
pected bias is now negligible at the cost of an increase in 
the statistical uncertainty. 

The dominant source of the increased statistical uncer- 
tainty is the remaining uncertainty in the amplitude of the 
patchy power spectrum. This uncertainty can be reduced 
by experiments sensitive to / > 3000 where the patchy 
power spectrum dominates, such as considered in the pre- 
vious Section. 
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Table 1 

Bias and standard deviation in the parameter estimation for Planck (^ max = 3000). 







Q. b h 2 


n s 


n' s 


m v 


Y He 




e s 


In Pi, 


w x 


bias 


0.00012 


0.00025 


-0.030 


-0.022 


-1.20 


-0.10 


1.9 


-0.0006 


0.03 


0.74 


la 


0.0013 


0.00030 


0.013 


0.0050 


0.45 


0.021 


3.2 


0.0026 


0.048 


0.56 



The parameters are: Q. m b? (matter density), VL b h? (baryon density), m v (neutrino mass), Yfj F (helium mass fraction), z r (rcionization 
rcdshift), 8 3 (angular sound horizon), w x (dark energy pressure to density ratio) and the primordial potential power spectrum is fc 3 Pj>(fc) = 
fc3p 4) (fc / )(fc/fc / ) n =- 1 + n = ln ( fc /' : /) with k f =0.05Mpc- 1 . 



Table 2 

Bias and standard deviation in the parameter estimation for a sample variance dominated experiment 

(£ max = 2000) 







0,bh 2 


n s 


n s 


m„ 


Yh c 


Zr 


0s 


InP* 


Wx 


bias 


-0.00018 


0.00030 


-0.030 


-0.0210 


-1.14 


-0.105 


0.88 


-0.00006 


0.019 


0.40 


1(7 


0.00100 


0.00027 


0.011 


0.0049 


0.43 


0.018 


2.87 


0.0020 


0.044 


0.52 



Table 3 

Ratio of bias to statistical uncertainty. 







n b h 2 


n s 


n' s 


m v 


Yh. 




0s 


In Pi, 


w x 


A 


bias/ o n 


0.027 


0.027 


-0.060 


-0.083 


-0.030 


-0.077 


-0.008 


-0.031 


-0.008 


0.007 


0.011 




1.0 


1.4 


1.1 


1.0 


2.4 


1.0 


1.0 


1.0 


1.1 


1.4 





Values arc for a sample variance dominated experiment (< maJ( = 3000) assuming the patchy signal has a fixed shape. Also shown is the ratio 
of the expected la errors (cr n ) to the previous values without the extra parameter (a ). 



6. IMPACT ON LENSING POTENTIAL RECONSTRUCTION 

To understand the secondary confusion related to patchy 
reionization in lensing reconstruction with CMB data, we 
follow suggested approaches in the literature (Hu & Oka- 
moto 2002; Cooray & Kesden 2003) and modify the noise 
calculation related to lensing reconstruction to include the 
kSZ contribution. In the case of temperature, the variance 
of quadratic statistics, with 2 (1) as the filtered version of 
the squared temperature (see, for example, Cooray & Kes- 
den 2003), leads to a reconstructed lensing potential power 
spectrum of 



(e 2 (i)e 2 (i')> = {2-k) 2 5 d {\ + 1') [of* + ni 

where the noise is 

d 2 h [l-l^ + l- (1 



bij 



(2tt) 2 



ZU h °|l-li| 



(29) 



(30) 



Here, Cf is the power spectrum of unlensed primordial 
temperature fluctuations, while, following equation (24), 
C\ at accounts for all contributions to the temperature power 
spectrum including noise. Instead of point sources in- 
cluded in equation (24), here, we include foreground con- 
tributions related to secondary effects, mainly thermal (tSZ) 
and kinetic (kSZ) SZ contributions. 

In Figure 8, we summarize our results where we show the 
deflection angle power spectrum, Cf d = 1(1 + l)Cf^ and 
related noise contributions. In addition to temperature, 
we have also considered lensing reconstruction from polar- 
ization data and have concentrated on the quadratic com- 
bination related to E and B-mode maps. This combina- 
tion provides the best estimator for lensing reconstruction 



from CMB data if the signal-to-noise is sufficiently high 
(Hu & Okamoto 2002). In calculating these noise curves, 
we assume temperature and polarization maps with noise 
rms of 2 [iK in each beam-size pixel, and a beam size of 
3'. Such sensitivities may be achieved by future ground or 
space-based missions. 

The thermal SZ contribution boosts the noise by an or- 
der of magnitude from the case with noise alone. This 
is, however, not a significant source of worry since the SZ 
thermal contribution can be removed from thermal CMB 
maps based on the frequency spectrum of the SZ contri- 
bution in multifrequency data. The source of worry will 
then be the kSZ signal, including the contribution from 
patchy reionization. We show the level of expected noise 
following our calculation related to Figure 1. The noise 
is a factor of a few higher than the case of detector noise 
alone. 

Equation (29) assumes that the patchy power spectrum 
is perfectly known. Uncertainties in the amplitude and 
shape will further increase the noise in the lensing poten- 
tial reconstruction, though probably not by much since 
measurements at I ~ 5000 combined with modeling (such 
as described in our paper) can be used to constrain the 
amplitude and shape. 

While the lensing reconstruction from temperature data 
is degraded, this is not the case for polarization. As shown 
in Figure 8, the noise level including patchy reionization is 
the same as when only detector noise is included. Even for 
a perfect polarization map with no detector noise, the re- 
construction with polarization is not significantly affected 
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Fig. 8. — The deflection angle power spectrum (black thick solid 
line) and the noise variance for its reconstruction with different as- 
sumptions. In solid red, we show the noise related to reconstruc- 
tion with temperature maps while in blue we show the same using 
polarization data (based on the quadratic combination of E mode 
with B mode). The dotted line shows the reconstruction noise with 
temperature in the presence of SZ thermal effect. This secondary 
contribution can be removed in multi frequency data using the SZ 
spectrum. The resulting lensing reconstruction is then affected by 
the kSZ signal from patchy reionization (red dashed line). The blue 
dashed line shows the same noise in the presence of kSZ polariza- 
tion contribution. The latter is significantly smaller and leads to 
the same noise when the polarization from patchy reionization is 
ignored. 



since the level of foreground noise is significantly below 
the cosmic variance limit of the E-mode power and be- 
low the cosmic variance of lensed scalars in the B-mode 
power. Thus, polarization observations provide a confu- 
sion free way to extract lensing information and to study 
primordial anisotropy contributions. 

In addition to foreground noise variance, given by their 
power spectra, the foreground secondary effects can also 
be correlated with lensing potentials. While this leads to 
an additional noise contribution, we have found this to 
be significantly smaller than the Gaussian noise shown in 
Fig. 1, since the quadratic maps are prefiltered to maxi- 
mize lensing reconstruction such that other higher order 
correlations are reduced (Cooray & Kesden 2003). An- 
other secondary source of noise is the non-Gaussianity of 
foregrounds themselves. In the case of kSZ related to 
patchy reionization, we expect this contribution to be neg- 
ligible since at high redshifts clustering is linear and on 
large scales we will be averaging over the contributions 
from many patches. 

7. SUMMARY 

The kSZ power spectrum is important phenomenolo- 
gically in two ways: 1) as a contaminant of the primary/len- 
sing CMB power spectrum and 2) as an interesting probe 
of the epoch of reionization. We have shown that, if ne- 
glected, the patchy contribution can lead to statistically 
significant biases in the cosmological parameters estimated 
from Planck data. Higher-resolution observations (such as 



by ACT, APEX, SPT x ) could be even more strongly af- 
fected. We have proposed a simple solution, which is to 
phcnomcnologically parameterize the patchy power spec- 
trum. Fortunately for this purpose, as we have seen an- 
alytically, although the amplitude in the relevant 1 range 
is highly uncertain, the shapes are not. For the models 
we have considered, extending our cosmological parame- 
ter set to include one more parameter (the amplitude of 
the patchy power spectrum) eliminates any significant bias 
(even for an all-sky cosmic-variance limited experiment to 
^ = 3000). 

In addition to parameter estimation, high-resolution tem- 
perature and polarization anisotropies are affected by grav- 
itational lensing. CMB maps can be used to reconstruct 
the lensing potential, whose auto and cross power spectra 
(cross with CMB maps and cosmic shear lensing potential 
maps) are sensitive to the dark energy properties and neu- 
trino masses. We have also investigated the degradation of 
lensing potential reconstruction and conclude that lensing 
studies with temperature data alone will be affected by 
fluctuations related to patchy reionization while polariza- 
tion data will not be affected. Thus parameter estimation 
and lensing-potential reconstruction from polarization are 
much less affected by reionization than is the case for tem- 
perature anisotropies. 

The kSZ power spectrum can also give us some infor- 
mation about the reionization process. At I > 3000 and at 
frequencies near the null of the thermal SZ effect, the kSZ 
reionization power spectrum will be the dominant source 
of fluctuation power on the sky, after cleaning out point 
sources. SZ survey instruments such as ACT, SPT, APEX 
and CARMA 2 with the ability to clean out the thermal 
SZ signal based on its spectral dependence will be able 
to make highly accurate measurements of the kSZ power 
spectrum. We have calculated the expected power spec- 
trum errors for a nominal observation, including model- 
ing of residual point source contamination. If r becomes 
tightly constrained, this type of experiment should allow 
us to determine the mean effective bias and thus provide 
important information about the sources of reionization. 

In this paper we have only calculated power spectra. We 
expect that these two-point functions are the most impor- 
tant statistics, both for consideration of patchy reioniza- 
tion as a signal and as a contaminant (see Castro 2003 for 
a discussion of non-Gaussianity from homogeneous reion- 
ization). At I < 10 4 , we expect the signal to be approxi- 
mately Gaussian due to the large number of patches along 
each line of sight. Non-Gaussian properties may also be 
of interest and have been studied by Gnedin & Shandarin 
(2002), though at very small scales (£ ~ 10 5 ). 
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2 Combined Array for Research in Millimeter- wave Astronomy 
(http:/ /www. mmarray.org) 
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